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Q\ [ Higher than the lowest Landau level contributions to magnetization and spe- 

cific heat of superconductors are calculated using Ginzburg - Landau equations 
approach. Corrections to the excitation spectrum around solution of these equa- 
tions (treated perturbatively) are found. Due to symmetries of the problem 
leading to numerous cancellations the range of validity of the LLL approximation 
in mean field is much wider then a naive range and extends all the way down to 
H = — 7g ■ Moreover the contribution of higher Landau levels is significantly 
smaller compared to LLL than expected naively. We show that like the LLL 



part the lattice excitation spectrum at small quasimomenta is softer than that of 
usual acoustic phonons. This enhanses the effect of fluctuations. The mean field 
calculation extends to third order, while the fluctuation contribution due to HLL 
is to one loop. This complements the earlier calculation of the LLL part to two 
loop order. 



I. Introduction 

Ginzburg - Landau effective description of high T c superconductors has been 
remarkably successful in describing various thermodynamical and transport prop- 
erties. 1 However when fluctuations are of importance, even this effective descrip- 
tion becomes very complicated. Some progress can be achieved when certain 
additional assumptions are made. One of the often made additional assumption 
is that only the lowest Landau level (LLL) significantly contributes to physical 
quantities of interest. 2 ~ 5 ' 7 ~ 9 There is a debate however on how restrictive the LLL 
approximation actually is. Naively when H < ^ (see the dotted line on Fig. 
1), even within mean field approximation, one should consider higher Landau 
levels (HLL) mixing in the Abrikosov vortex lattice solution of the GL equations. 
When fluctuations are included one can argue using Hartree approximation 10 that 
the LLL range of validity is even smaller. However direct application of the LLL 
scaling to magnetization and specific heat on YBCO suggest that the range of 
applicability is much wider - all the way down to 1 — 3T. 11 ' 8,12 It is not clear why 
HLL do not contribute. 

In this paper we explicitly calculate the effects of HLL at low temperatures 
in the vortex solid or liquid phase and establish the realistic range where the 
LLL approximation is valid (see the heavy dashed line on Fig.l). We reanalyze 
the HLL corrections to mean field equations going to higher order then in 13 and 
find that the expansion converges for H C 2 > H > ^. Importantly within this 
radius of convergence the LLL contribution constitutes more than 95%. Then we 
calculate the HLL fluctuation effects to one loop order complementing the LLL 
calculation to two loops by one of us 9 (later referred to as I). 

Ginzburg parameter Gi characterizing importance of thermal fluctuations is 
much larger in high T c superconductors then in the low temperature ones. More- 
over in the presence of magnetic field the importance of fluctuations in high T c 
superconductors is further enhanced. Under these circumstances corrections to 
various physical quantities like magnetization or specific heat are not negligible 
even at low temperatures. It is quite straightforward to systematically account for 
the fluctuations effect on magnetization, specific heat or conductivity perturba- 
tively above the mean field transition line using Ginzburg - Landau description. 16 
However in the interesting region below this line it turned out to be extremely 
difficult to develop a quantitative theory. 

Within LLL in order to approach the region below the mean field transition 
line T < T m f(H) Thouless 2 proposed a perturbative approach around homoge- 
neous (liquid) state was in which all the "bubble" diagrams are resumed. The 
series provide accurate results at high temperatures, but for the LLL dimen- 

sionless temperature a^ = ( GiT ^i H2 ) ™ < —2 become inapplicable. 

Generally attempts to extend the theory to lower temperature by Pade extrapo- 
lation were not successful. 4 Alternative, more direct approach to low temperature 



fluctuations physics is to start from the mean field solution and then take into 
account perturbatively fluctuations around this inhomogeneous solution. Experi- 
mentally it is reasonable since, for example, specific heat at low temperatures is a 
smooth function and the fluctuations contribution experimentally is quite small. 
For some time this was in disagreement with theoretical expectations. 

Eilenberger calculated spectrum of harmonic excitations of the triangular vor- 
tex lattice (see eq.(28) below) 3 and noted that the gapless mode is softer then 
the usual Goldstone mode expected as a result of spontaneous breaking of trans- 
lational invariance. The inverse propagator for the "phase" excitations behaves 
as k 2 z + const{k\ + k*). The influence of this unexpected additional "softness" 
apparently goes beyond enhancement of the contribution of fluctuations at lead- 
ing order. It leads to disastrous infrared divergencies at higher orders rendering 
the perturbation theory around the vortex state doubtful. One therefore tends 
to think that nonperturbative effects are so important that such a perturbation 
theory should be abandoned. 6 However it was shown in I that a closer look at 
the diagrams reveals that in fact one encounters actually only logarithmic diver- 
gencies. This makes the divergencies similar to so called "spurious" divergencies 
in the theory of critical phenomena with broken continuous symmetry and they 
exactly cancel at each order provided we are calculating a symmetric quantity. 
Qualitatively physics of fluctuating D = 3 GL model in magnetic field turns out 
to be similar to that of spin systems in D = 2 possessing a continuous symmetry. 
In particular, although within perturbation theory in thermodynamic limit the 
ordered phase (solid) exists only at T = 0, at low temperatures liquid differs very 
little in most aspects from solid. One can effectively use properly modified per- 
turbation theory to quantitatively study various properties of the vortex liquid 
phase. This perturbative approach agrees very well with the direct Monte Carlo 
simulation of. 8 The question arises whether one can extend the well controlled 
perturbative calculation beyond the LLL. Sometimes a hope is expressed that the 
additional softness is an accidental artifact of LLL approximation. The present 
work explicitly shows that this is not so. It is a fundamental general phenomenon. 

The paper is organized as follows. Model is described and perturbative mean 
field solution developed in section II. The expansion parameter will be the dis- 
tance from the mean field critical line a^ = |(l — f — j[~)- Range of validity 
of the expansion and of the LLL approximation is discussed. Then in section III 
we derive the spectrum of excitations to leading order and to the next to leading 
order in ah- The free energy to one loop is calculated in section IV. Section V 
contains expressions for magnetization and specific heat and discussion of validity 
range of the fluctuation contributions calculation. Finally we summarize the re- 
sults in section VI. Details of the mean field calculation can be found in Appendix 
A, while details of the HLL spectrum calculation can be found in Appendix B. 



II. Model and the perturbative mean field solution 

A. Model 

Our starting point is the GL free energy: 

Here A = (—By, 0) describes a nonfluctuating constant magnetic field. For 
strongly type II superconductors (k ~ 100) far from iJ c i(this is the range of 
interest in this paper) magnetic field is homogeneous to a high degree due to su- 
perposition from many vortices. For simplicity we assume a = a(l — t),t = T/T c 
although this dependence can be easily modified to better describe the experi- 
mental coherence length. 

Throughout most of the paper will use the following units. Unit of length 
is £ = \/h 2 / (2m ab aT c ) and unit of magnetic field is H c2 , so that dimensionless 
magnetic field is b = B/H c2 . The dimensionless Boltzmann factor in these units 
is (the order parameter field is rescaled as ip 2 — > ^j-^ip 2 ): 

f = lj d3 4 m2 + 1 19 ^ 2 - l -^^ 2+ ^ |4 ' (2) 

The dimensionless coefficient is 

to = V2GinH. (3) 

where the Ginzburg number is defined by Gi = |( 327re ^^P 7 — ) 2 and 7 = m c /m a 5 
is an anisotropy parameter. This coefficient determines the strength of fluctua- 
tions, but is irrelevant as far as mean field solutions are concerned. 

B. Mean field solution by expansion in ah 

Now we turn to a perturbative solution of the Ginzburg-Landau equations 
near the mixed state - normal phase transition line. This has been done before 13 
to the second order, however the range of applicability and precision of the LLL 
approximation at large n has not been fully explored. The z direction dependence 
of the solutions is trivial and will not be mentioned until fluctuations will be 
discussed. The expansion parameter is 

1-t-b 

^ = ;: • (4) 



Rewriting the quadratic part in terms of operator (" Hamiltonian" ) 7i = |(— D 2 — 
b) whose spectrum starts from zero, one obtains the following free energy density 
over T 



L = L = }l 

T uj uj 



d 2 x(i>*Hil>-a h \i>\ 2 + ^\i>\ 4 



The equation of motion is therefore 



(5) 



H^ - a h ij + ij\ij\ 2 = 
This equation is solved perturbatively in ah by assuming 



(6) 



$ = (a A ) i / i [$ + o h $ 1 



(7) 



It is convenient to represent $o> ^i, ••• i n the basis of eigenfunctions of 7i, 7i(p n = 
nb(p n , normalized to unit "Cooper pairs density" < |y? n | 2 >= J cell d 2 x\(p n \ 2 ^ = 1, 
where "cell" is a primitive cell of the vortex lattice. Assuming hexagonal lattice 
symmetry one explicitly has; 



•p 



J-nl n n\a *-^ a 

l=—oo 



exp < i 



7il(l - 1) 2nVb 



Ix 



2 a \ 



where A= = \/-%- is the lattice spacing. 



To order zero 



H% = 



(8) 



and $o is proportional to the Abrikosov vortex lattice solution tp which is the 
expression eq.(6) for n = 0: 



$ = go¥- 



To order k, one expands 



$i = &¥> + X^;v 



(9) 



n=l 



Inserting into eq.(6) one obtains to order a\: 



H<&i = g p - g \g \ 2 tp\tp\ 2 
Taking the inner product with tp one finds that 



(10) 



1 



9o 



w 



9l,n 



gr, 



nbpV 2 ' 



;n) 



where the Abrikosov's constant is the following average over primitive cell: f3 
(3a =< M 4 >~ 116 . Inner product with tp n determines g™: 



(12) 



where (3 n =< \tp\ 2 ip n ip* > . To find g\ we need in addition also the order a h 
equation: 



H$ 2 = $ 1 -(g ) 2 (2$ 1 \<p\ 2 + <S>* 1 <p 2 ) 
Inner product with tp gives: 

3^ ((3 n ) 2 



2 ^ n6/3 5 /2 ' 

n=l 



(13) 



(14) 



C. Mean field result for free energy. Orders o? h and a\ 

The mean field expression for the free energy to order a 2 h is well known. In- 
serting the next correction eq.(7) into eq.(5) one obtains the free energy density: 



j i 



mf 



T 



1 

UJ 



a h _ a h v^ 

OR R3h 2^ 



(P 



n\2 



2(3 [3 3 b 



n=l 



II 



1 



-.43a? - .0078-^ 
b 



(15) 



It is interesting to note that (3 n ^ only when n = 6j, where j is an integer. 
This is due to hexagonal symmetry of the vortex lattice. 13 For n = 6j it decreases 
very fast with j: (3q = —.2787, (3\2 = .0249. Because of this the coefficient of the 
next to leading order is very small (additional factor of 6 in the denominator). We 
might preliminarily conclude therefore that the perturbation theory in ah works 
much better that might be naively anticipated (see dashed line on Fig.l) and can 
be used very far from transition line. If we demand that the correction is smaller 
then the main contribution the corresponding line on the phase diagram will be 
b = .015 • (1 — t). For example the LLL melting line corresponds to a^ ~ 1. This 
overly optimistic conclusion is however incorrect as calculation of the following 
term in Appendix A shows. 
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D. Range of applicability of the expansion. How precise 

is LLL? 

Now we discuss in what region of the parameter space the expansion outlined 
above can be applied. First of all note that all the contributions to $1 are 
proportional to |. This is a general feature: the actual expansion parameter is 
%. One can as whether the expansion is convergent and, if yes, what is its radius 
of convergence. Looking just at the leading correction and comparing it to the 
LLL one gets a very optimistic estimate. For this purpose we calculated higher 
orders coefficients in Appendix A. The results for the $2 are following: 



n = J_ 

" l nb 



1 00 
& -JiJ2 9 ^ 2 < MM > + < 0,0\i,n >) 



(16) 



and 



1 
92 = 'W^ti ~2VB^2 9l9i{< °' ° lh j>+2< * 0| *' ° >} (17) 

where < ii,i2\ji,J2 >=< ip* 11 ^* 12 ^ 1 ^ 2 > and g 1 - when i = is defined to be 
equal to gj t ip % = tp when i — 0. 

We already can see that #2 an d #2 are proportional to (?™ and in addition 
there is a factor of 1/n. Since, due to hexagonal lattice symmetry all the g™ 
,n 7^ 6j vanish, so do g% ■ We checked that there is no more small parameters, 
so we conclude that the leading order coefficient is much larger than first (factor 
6-5), but the second is only 6 times larger than the third. 

The correction to free energy is 

T mf 1 0.056 4 



T u 6 2 b 2 v ' 

Accidental smallness by factor 1/6 of the coefficients in the % expansion due to 
symmetry means that the range of validity of this expansion is roughly ah < Qb 
or H < -^p Moreover additional smallness of all the HLL corrections compared 
to the LLL means that they constitute just several percent of the correct re- 
sult inside the region of applicability. To illustrate this point we plot on the 
Fig. 2 the perturbatively calculated solution for h — .l,t — .5. One can see that 
although the leading LLL function has very thick vortices (Fig. 2a), the first 
nonzero correction makes them of order of the coherence length (Fig. 2b). Fol- 
lowing correction of the order (^-) makes it practically indistinguishable from 
the numerical solution. Amazingly the order parameter between the vortices ap- 
proaches its vacuum value. Paradoxically starting from the region close to H C 2 



the perturbation theory knows to correct the order parameter so that it looks 
very similar to the London approximation (valid only close to H c i) result of well 
separated vortices. 

We conclude therefore that the expansion in cth/b works in the mean field 
better that one can naively expect. In the next section we investigate whether 
the same is true for the fluctuation contribution. 

III. Fluctuations spectrum 

A. Fluctuations to leading order in a 

To find an excitation spectrum in harmonic approximation one expands free 
energy functional around the solution found in the previous section. Within the 
LLL approximation this has been done by. 3 We generalize it to the case of all 
the Landau levels when perturbations due to nonlinear term are included. The 
fluctuating order parameter field ip should be divided into a nonfluctuating (mean 
field) part and a small fluctuation 

ij)(x) = $(x) +%(x). (19) 

The energy eq.(5) is then expanded in x retaining only quadratic terms 



h 



d x 



x*n x - a h \ x \ 2 + 2|$| 2 |xl 2 + i(**Y + $V 2 ) 



(20) 



Field x can be expanded in a basis of quasimomentum k eigenfunctions: 



V\ 



2tt 27T 

r ~ n , J2 H n(y^ + k x 1) 



l=—oo 



exp < i 



TV 



1(1 -1) My/te-j^ _ v - 



-ifev^ + t-^f) 2 



In addition instead of complex field Xk we wm use two "real" fields O k and A k 
satisfying O™ = 0*_ n k Ak = A™ k : 



xW 



X (x) 



n=0 
» oo 

Jk. „_ n 



(21) 



n=0 



In terms of these fields representing "optical" and "acoustic" phonons eq.(20) 
takes a form: 



» oo 

h = I Yl 2 ( nh - a h){O n k O n _ k + A n k A n _ k ) - 2a h (O k O_ k + A k A 

oo 

+ 2 J2 * k Al k Ki? + OiA>_ k L* + OL^M? + 0\Ol k N^ 

i,j=0 



(22) 



where elements of the matrix 



#' =< 1*1 v_v^ - *>>$) + ^*vv^ - *v_v*4) > 



We expand f 2 in a^. The order a^ term is 

» oo 

lj2-a h (OlOl k + AlAl k ) + 

Jk - i=0 

oo 

°* E [4^i*^t J "(l) + OiAi fc Lj*(l) + OL^M^!) + 0*0^(1 



(24) 



M 



We will use the degenerate perturbation theory similar to one used in Quan- 
tum Mechanics to calculate the correction to the eigenvalues of matrix the LLL 

states ( A,0 states) to order a\. The matrix H = I ,, ,. lis analogous to 

" Hamiltonian" while („) is analogous to eigenvector. Explicitly the matrix ele- 
ments are: 



(25) 



K J M = \< M V_V^ + vWi) + \i^^.t + ^V?) > 

W) = ^ < M V_V^ " *Wi) + \^^.r ~ <f<p-it<p*i) > 

K J (i) = ^ < M V_V^ - ?V<4) + ^*VV^ - *>V_V*4) > 

We will use definitions 

# =< ^WV^TP > ( 26 ) 

ik =< ^^fit^-Tt > 

When an index is zero we drop it throughout the paper. For example when n — 0, 
we call /3£ as /^ and when fc = 0, (3% = (3 n etc..). 

Considering to order a^ above matrix element i,j = of if, L, M, iV is: 



iYl = | /9fc+ 2^ Re7fc (27) 

1 
Li =M fc = --Im7 fc 



In deriving this we have used a property (3k = (3- 



k- 



We now diagonalize the matrix to find the eigenstates of LLL (which is of 
order a^). Eigenvalues are 



e A = a h [ -1 + -(3 k - -^|7fc| ) (2£ 



2^ 1, 
eo = a h [ -1 + -jPfc + ^|7fc| 
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as was found originally by Eilenberger. 3 The "acoustic" branch is shown on 
Fig. 3a. The rotation transforming in to these eigenstates is: 



A k = cos — A k + sin — O k 

~ u k u k 

O k = - sin — A k + cos — O k 



(29) 



where ^ k = \^ k \ exp[i9 k }. Similar calculation for n th Landau level gives the spec- 
trum: 



e n A ,o = a fc ( -1 + | < MV?V2 > =F^| < OTVZtfU > I ) (30) 



B. Spectrum of fluctuations beyond leading order in au 

In this subsection we calculate the correction of eigenvalues of LLL to order 
a 2 h . The " Hamiltonian" H in addition to has the o^ part H\ given in eq.(25) also 
have the a\ part H 2 . As will be explained in the next section, we will need only 
correction to the LLL to the a\ order, not the HLL. Therefore we will need only 
the i,j = matrix element of H 2 : 



K 2 



No 



oo _. 

Y — 

^ nbB 2 

71=1 



OO 

E 



L 2 + M 2 



^nb(3 2 

1 



oo 

E 



^nb(3 2 

n=l 



^(3 n 2 (2(3 k - Re 7fe ) - 2/3 n (Re/3£ + Re(3 n _ k - Re 7 ^ 
-(3 2 n {2(5 k + Re 7fc ) - 2(3 n {Repl+ Re^ + Re 7^) 



lf3 2 n l mik + ^ n l mi n 



(31) 



Note that we do not show L 2 and M 2 separately as our result will depend only 
on L 2 + M 2 . According to the degenerate perturbation theory we need to diago- 
nalize Hi which already has been done in the previous subsection and then use 
the resulting states A k and O k to calculate the second order correction to the 
eigenvalue:^ = a\{E,nag + E ff diag ). The diagonal contribution is 



E, 



diag 



K A k \H 2 \A k > 

= (cos h 2 K 2 + N 2 (sin ^) 2 + (L 2 + M 2 ) sin °-± cos 



&k_ 

2 



(32) 



Substituting the matrix elements eq.(31) we obtain 
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Pn r <jPn 

2[Re^ + Re^-cos^Re^-sin^Im^] } (33) 



E diag = J2 ^Up{^f( 2f3k ~ l 7fc D 

n=l 



In the off diagonal contribution , 

< AAHAn >< n\HAA k > 



J^offdiag ' / J 



nb 

n=l 



1^ n h^ 



nb ^os |) 2 [| < A^HMk > I 2 + I < MHi\O n k > | 2 ] (34) 



+ (sin|) 2 [| < O fc |ifi|Oj? > | 2 + | < O fc |ffi|A2 > | 2 ] + sin | cos | 
x[ < A fc |fli|^fc >< 4J|#i|O fc > + < A k \H x \O n k >< 02\Hi\O k > +c.c]} 

Details of calculation of these matrix elements can be found in Appendix B 
together with definitions of quantities F. The result is: 



Eoff d ia 9 = 4£^0**(1)| 2 + l^(l)| 2 + \F£(2)\ 2 + \FH k (2)\ 2 



+ -^m\l)\ 2 + |F_\(1)| 2 - |F fc "(2)| 2 - |F_\(2)| 2 ] (35) 



+ ^± lm[F n k{l)F1 n {2) _ Fr(l)F fc "(2)]} 

Similar result for O^can be obtained from the above formula by changing the sign 
of cos^and the sign of sin#fcin the formula above. 

It is crucial to see whether there is an k 2 term in higher orders for the " acous- 
tic" branch A. We calculated numerically the contributions to the spectrum till 
n — 8. All the k 2 contributions to any of them cancel. Moreover even all the k 4 
contributions for odd n cancel although the even n give negative contribution to 
the rotationally symmetric combination (k 2 + k 2 ) 2 . Numerically the coefficients 
are: 2.2 • 10" 6 , 5.0 • 10~ 5 , -6.3 • 1(T 6 , 4.7 • 1(T 7 for n = 2, 4, 6, 8 correspondingly. 
The resulting correction to the spectrum of " acoustic" branch due to the n = 2 
level is shown on Fig. 3b. 

After we have established the spectrum of elementary excitations of the Abrikosov| 
lattice, we are ready to calculate the fluctuations contributions to various physical 
quantities. 
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IV. Fluctuation contributions to free energy, magnetiza- 
tion and specific heat 

A. Higher Landau levels contributions to free energy 

The thermal fluctuation part is 



-T\og[z\ = r mf + r j 



fluci 



T 



Ffluc = ^J2\ Trln 



n=0 



e\{k) + | 



+ Tr In 



e n oik) + | 



(36) 



TL 2 x L z abJ2 



n=0 



< Je\{k) > + < Jel{k) > 



where a = i and we performed integration over k z . 

The LLL contribution to order y/a~h in 2D has been calculated by Eilenberger. 3 
The 3D result is 9 : 



•r(l/2) 
''flue 



aba 



1/2 



.(1) 



<V £ iW> + <VW> 



-(1)/ 



3.16(T ba 



1/2 



(37) 



We calculated it's higher ah correction which is of order a h ' using eq.(33) and 
eq.(35) 



T 



(3/2) 



flue 



-aba, 
2 h 



3/2 



.(2), 



-(2), 






-AAhaba 



3/2 



(38) 



As noted below e A (k) and e A (k) given in the last section contain contributions 
from mixing with all the HLL. Table 1 details contributions to this term from 
levels till n = 8. The contributions are negative for all n ^ 6j where j is an 
integer and positive otherwise. 

Table 1. 

Contributions to free energy of mixing of the LLL with HLL. 
Given in units of ^crba h ' . 



level n 


1 


2 


3 


4 


5 


6 


7 


8 


A mode 


-.253 


-.082 


-.053 


-.063 


-.063 


.247 


-.017 


-.005 


mode 


-.230 


-.086 


-.051 


-.023 


-.012 


.018 


-.003 


-.001 
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The contribution of HLL is: 

oo 



•J flue 

T 



n=l 

oc 



, Jn 

n=l n=\ 



< Jnb + a h e n A {k) > + < Jnb + a h e^(k) > 

1 °° 1 

2ab^ 2 J2^+ ^ 2 a h J2^i< £ A( k ) > + < £ o(k) >] + 0{a\) (39) 



- IfT(O) , -r-(l) \ 

— rp W flue "T ''flue) 

The first divergent (as powers 3/2 and 1/2 in the ultraviolet) term renormalizes 
energy. However it has a finite magnetic field dependant part which should be 
calculated by subtracting the 6 = value of the free energy. The proper regu- 
larization is made by restricting the number of Landau levels and then showing 
that after regularization the answer does not depend on it. The calculation is the 
same as one done in the normal phase (obviously UV divergencies are insensitive 
to the phase in which they are calculated), see for details and discussion. 14 The 
result is: 



T 



(0) 



flue 



T 



.526a6 3/2 (40) 



This exhibits the diamagnetic nature of the bosonic field. 15 The second term is 
proportional to ^ A= and also diverges but only as power 1/2 in the ultraviolet 
and renormalizes a^. To see this we calculate the sum 

< e\{k) + e n (k) > k = -2 + i « M VITVZ >*>*= -2 + j 

The last equality follows from the curious property of (p* k n (x)(p^(x) that it depends 
only on XiVb — £ijkj. We see that apart of renormalizations there is a finite 
correction: 

J flue - -- / - *• 



1.459 (-1 + - j ab 1/2 a h (41) 



T 
The following one is of the order a? h and will not be calculated here. 

V. Results for magnetization and specific heat. Range of 
applicability of the loop expansion 

Here we discuss the nature and range of applicability of the expansions we 
used for fluctuating superconductors (for which Gi is not negligibly small). There 
are two small parameters used. The first one is ^ which controls the expansion 
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of the mean field solution and therefore the HLL corrections was already discussed 
in section IID. The second small parameter controls the fluctuations. We assumed 
the mean field is the leading order and then expanded the statistical sum around 
it. Summarizing all the corrections the free energy density is 



T_ 
T 



= uj- 1 


2 h 


c 2 +c 3 T + c 4 ^ — j 


+ 


aba]l 2 


" (o) (a h \-y* jo) (o) /a h \ V2 (0) a h ~ 

c ° Vb) +c v2 + c i [t) +c 3/2y 


2 1,2 

ua o a 


-1 

h 


c (1) " 







+ 



(42) 



(43) 



where the coefficients (upper index is the power of u and the lower index is the 
power of ah) are: 



c (_1) 

c 2 

c (0) 

C l/2 



-0.434,4 X) 



-0.0078, c£ 0) = 0.526, 



(44) 



316, cj 0) = 1.06,4/2 



-0.445, cS = -118 



and a> is defined in eq.(3). The last term is the two loop contribution calculated 
in I. One clearly see that ooba h always appears together with an important 
"loop factor" a = 3 j 2 ~ 0.11. The expansion parameter therefore is 



auoba 



-3/2 



7rV2Gi 



tb 



1 



(1-t-b) 



3/2 



^larl 3 / 2 



In the last equation a^ is the often used dimensionless LLL temperature intro- 
duced by Thouless. 2 For Gi = 0.01 the condition cruoa^ ' < 1 is represented by 
the area above the dotted line on Fig.l 

Correspondingly the scaled magnetization is m = — ^£: 



m 
T 



1 



1 I J" 1 ) 



.(0) 



a h + 



^r { - l) h- l n 2 4- 9r> (_1) /,- 2 n 3 ^ 4 n"^ (_1) , r^n 3 3 r^h 1 / 2 r (0) n 1/2 

-c 3 % + 2c 4 o ft + o\c 3 cj a h - -c < - c 1/2 a h 



+ ^ba- 1/2 - \cTb-^a h + \c^ + 3 c (o) a i /2) _ ^ f^ + 1^., 



c« 



(45) 



»2 

while the scaled specific heat is c = —t^T: 



- = tT c u) 
t 



-i 



(- 1 ), ,~1 



-—Co o; 



2 

2 C 3/2 L 



J. r {0) /) 1 /2 , 3 J°) 1/2 3 _(0) /n -l/2\ , T 2r : 

+ C l ° + o C 3/2 a /i ~~ ^6 3/2 h 



34 1} & V - 94 1] b 2 a 2 h \ + aT c (-c ( ^ 2 ba h 1/2 + ^cfj 2 bta' h 

(46) 



-a». 
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VI. Conclusion. 

In this paper we showed why the LLL results are often valid far beyond the 
naive limit of applicability of the approximation for both the mean field and the 
fluctuation parts. Our results are valid strictly speaking between long dashed 
line representing H = c ^ ; and one of the dashed curves indicating the range of 
validity of the loop expansion for the fluctuation contribution (depends on value 
of the Ginzburg number Gi) . For nonfluctuating strongly type II superconductors 
our results can be directly checked by experiments done at low temperature or 
numerical solution (or even the "London limit approximation") and are in clear 
agreement. For small, but not very small Ginzburg parameter Gi one can compare 
with existing Monte Carlo simulations 8 ' 17 or experiments. Of course one can use 
existing high temperature expansion 2 to interpolate to the present expansion 
range. Results for LLL were presented in I and the HLL do not alter them 
significantly. The agreement with the MC simulations is very good although 
obviously the melting transition is not seen. As argued in I it is not expected 
to exist within the present model. The HLL do not change this conclusion. On 
the contrary we explicitly showed that the " supersoft" A mode has a propagator 
l/(k% + const{k\ + ky)) beyond the LLL approximation laying to rest a suspicion 
that this is a fluke due to LLL. This indicates that this unusual "softness" is due 
to some underlying symmetry which have yet to be explicitly identified. 
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Appendix A: Third order correction to the mean 
field solution and free energy 

In this appendix we provide some details of the third order in ah calculation 
of the mean field solution of the GL equations. 

To calculate (?2,one takes the inner product of Lp n on the two sides of eq.(13) 
and obtains eq.(16). To calculate #2, we need to consider the GL equation to 

7/2. 

ft$ 3 = $2 - [($o) 2 $; + ($i) 2 $o + 2|$o| 2 $2 + 2|$i| 2 $5] (47) 
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order a h 



Scalar product with ip gives eq.(17). 

Now we compute the a\ order correction to free energy. Substituting $2 from 
eqs.(16) and (17) we find: 



-a\ (< $ 2 |^|$i > + < $i|7i|$ 2 > - < $21^0 > - < $i|$2 > - < $i|$i >) 



(4£ 



92 -a 3 
> a h 



/3V2 



+«*E 



n=0 L 



nbtfti-zW 



n\2 



Dominant contributions come from: < 6, 6|0, >=< 0,0|6,6 >= .80260 
< 6, 0|6, >= .80283and those coefficients are real. 



Appendix B: Second order correction to 
fluctuations spectrum 

In this appendix we list matrix elements of the correction Hi given by eq.(25) 
between various states used in the calculation of the second order correction to 
energies of excitations. 






1 



on 1 o*n 1 ( -,n i _,*tA 

P-k + Pk + g l^-fe + 7k J 



(49) 



< Akl^iio^ > = - < M 2 ( w f - <P n _Tt<plit) > + 



2(3 



< <Aw% - SvUvr > 



1 
1 



on 1 o*n 1 f-,n -,*n\ 

-p_ k + p k + - (7-fe - 7 fc ) 



(50) 
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<O k \H 1 \A^> = - 
<Al\H 1 \O k >= 1 - 
< A k \Ih\Ol > = % - 
<Ol\H l \O k > = ^ 



on o*n i I 'n „,*n\ 

P- k - p k + - (7-fe - 7 fc ) 

on o*n , / „,*ra , „,n\ 

p k -p_ fc + - (-7-fe + 7fcj 
P k - P_ fc + - (-7-fc + 7fc j 



o*n i on i f^,*n , _,n\ 

P- k + P k + o (7-fc + 7 fc j 



(51) 



etc. From those formula, we can show 

|<0J^|0™>i 2 + |<0Ji?;u">| 2 



ff 2 

2 

J 2 



F k (l)\ 2 + \F- k (l)\ 2 ] 
[l^(2)| 2 + |F_\(2)| 2 ] (52) 



< A k \HAA n k >< A n AHAO k > + < A k \HM >< Ol\HAO k > +c.c 



2 



- 2 i F ra)Fm - j?*(i«(2)] + c . c , 



where 






(53) 



Finally we can show that 



E, 



offdiag 



■\ E ^ DW)! 2 + l^(!)| 2 + l^(2)| 2 + l^-\(2)| 2 ] 



^[l^(l)| 2 



+ 



sin#fc 



i^- fc (i; 



l^"(2)P 



-*(2)| 



2Im[F^(l)F™(2)-Fr(l)F fc "(2)]} 



(54) 
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Figure captions 



Fig. 1 



The range of the validity of the expansions in ah and the loop expansion. The 
region above the dotted line is the naively expected validity range of the LLL 
approximation. The region above the long dashed line is the actual validity range 
for the expansion of the mean field equations. The loop expansion applicability 
range lies below the dashed curves. We plot two curves with different values of 
Ginzburg number Gi, Gi = 0.1 and Gi = 0.01. The validity combining the mean 
field expansion and the loop expansion lies therefore between the long dashed 
line and the dashed curves. 



Fig.2 



The density |-?/> 2 |of the mean field Abrikosov solution for b = 0.1, t = 0.5. 
Fig. 2a is the lowest order approximation (LLL). Fig. 2b is the solution with the 
next order correction included, while in Fig. 2c the next next order correction is 
included. 



Fig.3 



The shear mode A spectrum. Fig. 3a is the spectrum obtained within the LLL 
approximation. Fig. 3b is the correction to the spectrum when the HLL mixing 
effect is considered. 
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